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Abstract 

We present data on pp tjtjtt^ at beam momenta of 600, 900, 1050, and 1200 MeV/c. 
At the higher three momenta, a signal is clearly visible due to pp /o(1770)7r'^, /o(1770) — >■ 
ijr]. It has mass 1770 it 12 MeV and width 220 it 40 MeV, where errors cover systematic 
uncertainties as well as statistics. 

PACS: 13.75. Cs; 14.20. GK; 14.40 
Keywords: Mesons; Resonances; Annihilation 

We investigate pp -> rjrjiT^ at p momenta from 600 to 1200 MeV/c. Data were taken at 
LEAR by the Crystal Barrel Collaboration. Data from 1350 to 1940 MeV/c are reported in 
an accompanying paper p. They provide evidence for /o(2100) rjrj, a2(1660) — > rjn and a 
broad /2(1980) ^ 7777. The first resonance lies above the mass range of present data, but is an 
important element in 0^ spectroscopy. The 02(1660) is strong at higher momenta but, because 
its threshold is at ~ 2200 MeV, it makes only a minor contribution to present data up to 1200 
MeV/c. The /2(1980) plays an important role, since it is broad and extends into the present 
mass range. We shall provide evidence for an / = 0, J^'-'O^'^ resonance at 1770 MeV, decaying 
to rjrj. This resonance is important in trying to unravel the mysteries concerning the /j(1710) 
[2], hence in establishing the systematics of 0"*" mesons. 

At lower masses, there has been earlier evidence from the Crystal Barrel collaboration for 
an / = resonance whose mass is variously estimated from 1300 to 1380 [3-6]. There is also 
an / = 1 resonance at 1450 MeV [7] and finally the well known /o(1500). At higher mass, the 
WA102 collaboration has published evidence for an /o(2020) resonance decaying to 47r [8]. 



Momentum 


Events 


Background 


(MeV/c) 




(%) 


600 


2922 


2.8 


900 


9023 


3.9 


1050 


6607 


3.3 


1200 


9594 


3.7 



Table 1: Numbers of events and estimated background levels 
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The experimental conditions and data processing are identical to those of the higher momenta 
reported in Ref . pQ , and we refer to that paper for experimental details. The largest backgrounds 
are from r77r°7r°7r° and cutt^'tt", two photons being lost in the former case and one in the latter. 
Numbers of events and estimated backgrounds are shown in Table 1. Using cross sections 
observed in these other channels, we find that backgrounds are close to a phase space distribution. 
They are included in fitting the data, but omitting them has negligible effect on the analysis, 
where log likelihood changes by < 4. 

Fig. 1 shows Dalitz plots at the four beam momenta and mass projections on to rji] and 
rjiT. Histograms superimposed on the projections are the result of the maximum likelohood 
fit described below. There are conspicuous peaks due to /o(1500), ao(980) and a2(1320). The 
feature of interest is a shoulder at ~ 1770 MeV in rjr], clearly visible in Figs. 1(f), (g) and (h) 
at 900, 1050 and 1200 MeV/c. It is barely visible at the higher momenta reported in Ref. [1]. 
At 1350 MeV/c, there is a shoulder at ~ 1770 MeV, and the amplitude analysis gives a small 
optimum there, but the significance is found to be much less than at lower momenta. 

We shall report here a full partial wave analysis of production and decays of resonances. 
Amplitudes we find to be significant are listed in Table 2. In making this choice, we are guided 
by experience that it requires roughly 225 MeV/c of momentum per unit of orbital angular 
momentum L in the production process. Present data are consistent with this guideline. The 
last column of Table 2 shows momenta in the decay process at the highest beam momentum of 
1200 MeV/c. For the ao(980)77 channel, we require up to L = 3. There is no evidence for the 
presence of L = 4; this agrees with the expectation that the first 4~ qq state will lie at a mass 
of 2300 MeV or above, i.e. a beam momentum of at least 1650 MeV/c. The ao(980) is fitted 
with a Flatte form using parameters derived from earlier data on pp — > rjTT^Tc^ at rest P] . The 
ao(980) signal is slightly stronger in the data than the fit, particularly in Fig.l(i); the fit can be 
improved by reducing the width of the ao(980) by ~ 20%. The /o(1500)7r channel is described 
successfully with waves up to L = 2. The strong a2(1320)r7 channel requires up to L = 3, but 
again there is no evidence for L = 4; the 2"*" initial state with helicity makes quite a significant 
L = 3 contribution: (4 — 8)% at 900 MeV/c and above. 

Channel Partial waves CM momentum 



(MeV/c) 



ao(980)r/ 




734 


/o(1500)7r 


^'S'o,'^Pi,^-D2 


539 


02(1320)7/ 


^2(1 = and 2),^So{L = 2) 
3Pi, 2+(m = and 1), ^F^iL = 1) 


498 








2+(L = 3,m = 0) 




/o(1770)7r 


^5*0, 3Pi 


323 


/2(1270)7r 




690 


/2(1980)7r 


^D2?Pi,2+{m = and 1),^F3 


100 



Table 2: Partial waves included in the analysis; the third column shows the centre of mass 
momentum on resonance for the highest beam momentum of 1200 MeV/c 

Figs. 2 and 3 show examples of angular distributions in bands centred on resonances. The 
fits, shown by histograms, are equally good at other momenta. Fig. 2 shows the production 
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Figure 1: (a)-(d) Dalitz plots at 600, 900, 1050 and 1200 MeV/c, (e)-(h) Projections on to 
M{rjri), (i)-(l) Projections on to M{ri7r). 
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Figure 2: A comparison of the fit at 1200 MeV/c with angular distributions (unnormahsed) 
against cos r, where r is the centre of mass angular production angle; (a) a 65 MeV wide band 
around ao(980), (b) a 118 MeV wide band around a2(1320), (c) an rji] band from 1640 to 1840 
MeV, (d) a 120 MeV wide band around /o(1500). Points with errors are data and the histogram 
shows the fit. 




angular distribution in terms of the centre of mass angle r at which the resonance is produced. 
Fig. 3 shows decay angular distributions of a resonance, e.g. 02(1320) — ?■ rjn, in terms of its 
decay angle a with respect to the beam direction. 

Our data on the 37r° state (not shown here) contain strong /2(1270)7r° signals and lead us 
to expect weak /2(1270) rjr] contributions at the level of 3 — 5% in present data; these are 
barely visible in the projections. We find that this weak channel can be approximated with 
L = production only; this approximation has no impact on the interesting rjr] mass range 
around 1770 MeV. We have tried including /o(1370) in the fit, but no significant contribution 
is required. 

We find that most partial waves follow within errors a smooth dependence on centre of mass 
momentum, but with fiuctuations which are typically 25% for the larger waves (~ 8% branching 
ratio) and 50% for smaller waves contributing branching ratios of only 2%. Contributions below 
1% are dropped. We have data on pp — >■ rjTr^TT^ with roughly a factor 12-16 higher statistics 
|10] . There, many resonances are visible, but phase variations with momentum are mostly small, 
i.e. resonances are phase coherent. In present data, phase variations are again small, although 
there are variations of typically ±(10 — 25)° from one momentum to another. The fit is made 
by the maximum likelihood method. Our definition is such that a change of log likelihood of 
0.5 corresponds to one standard deviation for a change of one parameter. The acceptance of 
the detector is simulated by Monte Carlo events, processed through exactly the same selection 
procedure as data; statistics of accepted Monte Carlo events are roughly 4 times those of data 
events. 




0,2 0.4 0.6 0.8 1 0.2 0.4 0.5 O.S 



fj(1700)cosa f0(1500) cos a 




o0(980) cos a g2(1 320) cos a 



Figure 3: As Fig. 2 for decay angular distibutions with respect to the beam in the rest frame 
of the resonance. 



The fitting algorithm is sufficiently fast that the amplitude analysis produces a fit in roughly 
60 seconds of computing. This is sufficiently fast that we have been able to examine a large 
number (> 1500) of trials with varying ingredients in order to study the systematics of the fit, its 
stability and sensitivity to every channel. Although some ambiguities arc present at individual 
momenta (particularly between L = and L = 2 and between L = 1 and L = 3), we have 
located only one solution having a smooth variation of magnitude, and particularly phase, with 
momentum. The Dalitz plots of Fig. 1 are reproduced so closely by the fit that no difference 
from data is visible within statistics. 

The rji] signal at 1770 MeV is visible in rjrj projections of Figs. 1(f), (g) and (h) at 900, 1050 
and 1200 McV/c. Fits without any resonance in this mass range are shown in Figs. 4(a)-(c). 
There is an obvious surplus of events between 1700 and 1800 MeV. We have tried fitting it 
with J = or 2. Variations of log likehhood with its mass are shown in Figs. 5 (a)- (c). There 
is a strong, well defined optimum for J = (full curve). Optimum masses are 1758 ± 10 MeV 
at a beam momentum of 900 MeV/c, 1770 ± 8 at 1050 MeV/c and 1775 ± 8 MeV at 1200 
MeV/c. Branching fractions are respectively 5.5%, 11.7% and 10.1% at these three momenta. 
At 1350 MeV/c, there is a shallow minimum at the same mass and width, but the fitted /o(1770) 
contribution is only 1.5%. Data at 600 MeV/c are too close to threshold to show an optimum 
in the mass or width of /o(1770). Column (a) of Table 3 lists changes in log likehhood when 
/o(1770) is omitted from the fit. 

Spin gives a better optimum than J = 2 (shown by dashed curves in Fig. 5). Difi'erences of 
log likelihood are given in column (b) of Table 3. However, spin 2 still shows some optimum at 
a similar mass to spin 0. We find this is due to feed-through from mis-identified spin signal. 
The detector has holes around the beam with an opening angle of 12°. This incomplete coverage 
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Figure 4: Fits without any /j(1770) at beam momenta of (a) 900, (b) 1050 and (c) 1200 MeV/c. 
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Momentum 


(a) 


(b) 


(c) 


(d) 


(MeV/c) 










900 


-81.4 


-33.8 


16.0 


3.1 


1050 


-101.6 


-34.2 


19.3 


4.8 


1200 


-83.0 


-28.9 


17.5 


3.5 


1350 


-7.7 


-5.4 


19.2 


3.4 



Table 3: Changes to log likelihood of the basic fit, (a) dropping /o(1770), (b) replacing /o(1770) 
by /2(1770), (c) adding /2(1712) or (d) /o(1712) 



allows some cross-talk between spin and spin 2. In order to examine this, we have adopted two 
simulation procedures. In the first, we have used Monte Carlo events to generate data samples 
with the same partial wave amplitudes as are fitted to /o(1770) plus other channels. Then we 
have fitted these samples replacing /o(1770) by a 2+ resonance. This simulation reproduces the 
magnitudes of the dashed curves of Fig. 5(a) and (b) at 900 and 1050 MeV/c quite well. 

The fits with spin 2 use five partial waves ^1^2, '^-Pi, 2~^{m = and 1) and ^F^, while fits 
with spin use two partial waves ^5*0 and ^Pi. Statistically, we would expect the fits with 
spin 2 to be better than those with spin by 3 in log likelihood, because of the 6 extra fitting 
parameters. In practice, we find this to be an under-estimate, since certain combinations of 
amplitudes producing a J = 2 resonance with L = and 1 can simulate closely the flat angular 
distribution of the dominant ^5*0 amplitude producing J = 0. To examine this, our second 
simulation is to fit fake J = 2 resonances of width 220 MeV to the simulated data sets. What 
we observe in 10 such trials is that the fake 2"*" signal improves log likelihood by an average of 
17 ± 3. This sets a scale of what may be fitted incorrectly by spin 2; any improvement smaller 
than this cannot be believed. Our conclusion is that a fake 2"*" signal produced with L = or 1 
from 5 initial partial waves can be expected to fit the data better than a fake J = 0, produced 
from 2 initial partial waves, by approximately 12 in log likelihood. If we use this result, the 
conclusion is that J = is preferred over J = 2 by ~ 46 in log likelihood at 900 MeV/c, by 
~ 46 at 1050 MeV/c and by ~ 41 at 1200 MeV/c. These are statistically each over 9 standard 
deviations. If, pessimistically, we ignore the difference in the number of fitted partial waves and 
use just raw differences in log likelihood between spin and 2, the discrimination is still over 8 
standard deviations at each momentum. 

Fig. 5(d) shows log likelihood versus the width of the /o(1770) at 1050 and 1200 MeV/c, 
where the data are most definitive. At other beam momenta, the optimum is somewhat less 
well defined and flattens out at large widths. From our studies with a variety of ingredients, 
we estimate that the systematic errors on the mass and width of the resonance are respectively 
±12 and ±40 MeV. 

There has been previous evidence for /o(1770). A small peak has been observed at 1750 
MeV for J = in two sets of data on vrvr KsKs [11,12]. Secondly, a partial wave analysis 
of J/^ -> 7(47r) has fitted 47r peaks at 1780 MeV with J = [13]. The GAMS collaboration 
has reported [14] an X(1740) decaying to rjrj at 1744 ± 15 MeV, but they find a width < 90 
MeV. Their angular distribution is fiat, but spin is not claimed because of uncertainty about 
the production mechanism. E760 data [15] on pp — )■ rirjn^ show peaks at 1500, 1748 ± 10 and 
2100 MeV; they do not report a spin-parity analysis. In view of the fact that we observe strong 
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Figure 5: Variations of log likelihood with the mass of /o(1770), full line J = 0, dashed line 
J = 2, (a) at 900 MeV/c, (b) at 1050 MeV/c, (c) at 1200 MeV/c' (d) the variation of log 
likehhood with width at 1050 MeV/c (full curve) and 1200 MeV/c (dashed). 



production of /o(1500) and /o(2100) in our r/r/vr'^ data, it seems plausible that their peak at 
~ 1750 MeV may be the same resonance as we report here. Their width of 264 ± 25 MeV is 
acceptably close to ours. The BES group observes an /q at 1781 MeV in J /ip — )■ '-^{K^K~) |T6]. 

Two 0"*" qq resonances are to be expected in this mass range. One is the radial excitation 
of /o(1300 — 1380). The second is the missing ss state from a 0+ nonet made up of /o(1300 — 
1380), ao(1450) and i^*(1430). The resonance we observe could be either of these or a linear 
combination. The L3 collaboration [T7] has reported a narrow signal at 1793 ±18 MeV in KsKs, 
but without a spin-parity determination. A more recent analysis gives M = 1770 ±20 MeV, 
r = 235 ± 47 MeV, very close to values reported here. 

There is a long-standing debate over the spin of /j(1710), reported by the Particle Data 
Group as having a mass of 1712 ± 8 MeV and a width of 133 ± 14 MeV. The mass we observe, 
1770 MeV, seems too far removed from 1712 MeV to be the same resonance. It seems likely that 
there is one resonance at 1710 MeV with J = or 2 and a second at 1770 MeV. This is what 
the BES collaboration found [16]. There is then the possibility that the /j(1710) is also weakly 
present in our data. At 900 MeV/c, there is a small surplus of events in the 7777 projection of 
Fig. 1(f) from 1600 to 1700 MeV. At 1200 MeV/c there is again some indication of extra events 
in four bins close to 1700 MeV. Our fits include the 500 MeV wide /2(1980) demanded by the 
data of Ref. [1] at higher beam momenta. In ref. [1], the 2"*" signal is visible over the whole mass 
range from 1550 to 2000 MeV. The cos a dependence of Fig. 3(a) at 1200 MeV/c is again clear 
evidence for the presence of non-zero spin, but it does not appear to be narrowly concentrated 
in mass. Again, it is to be associated with /2(1980). At lower beam momenta, where the phase 
space is limited, there is the possibility of confusion between /2(1980) and /2(1710). At all beam 
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momenta, wc have tried adding to the present fit /o(1712) or /2(1712) with the PDG width of 
133 MeV. Improvements to the fit are hsted in columns (c) and (d) of Table 3 respectively. They 
are not significant, compared with our test with fake signals; they make no visible improvement 
to mass projections. If one tries to fit the extra events of Fig. 1(h) with /j(1710), the optimum 
mass is 1690 MeV but the width required is < 40 MeV and optimises at 25 MeV. This width 
is so much less than the PDG value for /j(1710) that we have no confidence in the presence 
of an additional /j(1710). An upper limit (90% confidence level) on its branching ratio at all 
momenta is 4%. 

In summary, there is strong (> 8a) evidence at each of three momenta for an /o with M — 
1770 ± 12 MeV and F = 220 ± 40 MeV. There is no significant evidence for a further resonance 
at 1710 MeV. 
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